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We discuss the distribution of ions around highly charged PEs when there is competition between 
monovalent and multivalent ions, pointing out that in this case the number of condensed ions is 
sensitive to short-range interactions, salt, and model-dependent approximations. This sensitivity 
is discussed in the context of recent experiments on DNA aggregation, induced by multivalent 
counterions such as spermine and spermidine. 
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INTRODUCTION 



Despite extensive theoretical and experimental re- 
search, polyelectrolyte (PE) solutions are relatively 
poorly understood compared to their neutral counter- 
parts Jj. The main difficulty in their theoretical treat- 
ment arises from the long-range nature of electrostatic 
interactions between the charged groups along the PEs. 
Another major difficulty arises in highly charged PEs 
due to their coupling with the surrounding ionic solu- 
tion, which is difficult to treat theoretically, since one 
cannot simply trace over the ionic degrees of freedom via 
the linearized Debye-Hiickel theory. 



In this paper we address the distribution of ions near 
highly charged PEs, concentrating on the case where 
more than one counterion species is present in the so- 
lution. We point out that the number of condensed ions 
is then highly sensitive to various parameters such as 
short range interactions, salt concentration, and model- 
dependent approximations - even at low concentrations 
of salt. In contrast, in solutions with only one type of 
counterion these parameters are important at high salt 
concentrations, e.g., in the ion-dependent solubility of 
proteins; at lower salt concentrations, typically up to 100 
mM, their influence on ion condensation is weak. 



After illustrating the above points using a simple ex- 
ample (Sec. 2), we discuss the competition between 
monovalent and multivalent ions in the context of DNA 
aggregation (Sec. 3), concentrating in particular on the 
role played by short-range interactions in the dilute (non- 
aggregated) phase. In Sec. 4 we discuss qualitatively the 
dependence on salt concentration of the number of con- 
densed ions. For this purpose we use a simplified two- 
phase model similar to Manning's model. 



2. CONDENSATION ON A SINGLE 
POLYELECTROLYTE CHAIN 

Let us consider first the distribution of ions near a 
single PE, i.e., taking the limit of infinite dilution for the 
PE solution. Let us assume also that the PE is uniformly 
charged with a charge per unit length equal to p. Suppose 
first that there is only one type of counterions in the 
system. 

When there is no salt in the solution only some of the 
ions remain bound to the PE, while the others escape to 
infinity. The number of bound ions, per unit length, is 
given by the well-known formula obtained from Manning 
condensation theory 2]: 
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where p* — 1/(zIb) and Ib — e'^ /{skBT) is the Bjerrum 
length, z is the counterion valency, e is the unit charge, 
e is the dielectric constant of the solvent and fesT is the 
thermal energy. When p > p* the condensed ions par- 
tially neutralize the PE such that its effective charge per 
unit length is equal to p* . This result is a consequence 
of the interplay between entropy and electrostatic energy 
at large ion-PE distances. Hence short range interactions 
(ion-ion or ion-PE) are not expected to modify the num- 
ber of condensed ions, although they may influence the 
distribution of the ions within the condensed layer. 

The number of ions in the vicinity of the polymer is also 
insensitive to salt at low and moderate concentrations |3] , 
as illustrated in Fig. 1 a. The figure shows the accumu- 
lated number of multivalent ions (^=4) per unit length up 
to a distance r from a PE having roughly the parameters 
of DNA (charge per unit length equal to 1/(1.7 A) and a 
radius of 10 A), as a function of r. The ion distribution 
is modeled using mean field theory and calculated using 
the Poisson-Boltzmann (PB) equation. This is done for 
simplicity, while in fact correlation effects beyond mean 
field are important in this case since the ions are multiva- 
lent. The dashed line shows the distribution when there 
is no salt in the solution; as r increases the accumulated 
number of ions Pb(r) approaches a constant, which is the 
number of condensed ions per unit length predicted in 
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FIG. 1: Accumulated number of 4-valent counterions, pbif), up to a distance r from DNA, modeled as a uniformly charged 
cylinder of radius 10 A: In (a), with only 4:1 salt of concentrations 0.01, 0.1, 1, 10, and 100 mM (solid lines). The limiting case 
of no salt (only 4-valent counterions) is also shown (dashed line). In (b) 10 mM of 1:1 salt is also present in the solution (solid 
lines). Dotted lines correspond to larger 4-valent counterions, having a distance of closest approach of 12 A to the DNA (the 
arrows connect solid and dotted lines corresponding to the same 4:1 salt concentrations of 0.01, 0.1, and ImM). The dashed 
line is the same as in part (a). The right axis in both figures shows azpb{r), the part of DNA charge compensated by the 
multivalent ions. In both parts the ion density profiles are calculated numerically using the Poisson-Boltzmann (PB) equation. 



Eq. 1^. When 4:1 salt is added to the solution at in- 
creasing concentrations oi Cz = 0.01, 0.1 1 and 10 mM 
(solid lines) the number of counterions close to the poly- 
mer is almost unaffected. A significant effect is seen only 
with Cz = 100 mM. 

Short-range interactions have no effect on the num- 
ber of condensed ions in the limit of zero salt concen- 
trations, and their effect remains small at low salt con- 
centrations. For example, the effect of ion-PE disper- 
sion forces was recently estimated for condensation of 
monovalent ions on DNA, using the Poisson-Boltzmann 
equation and the Hamaker approximation for the ion-PE 
dispersion interaction j^. The effect was considerable at 
1 M salt concentration, where electrostatic interactions 
are highly screened. However below 100 mM, even with 
rather strong dispersion interactions, in the order of four 
times the thermal energy close to contact, the effect on 
ion condensation was small. 

Returning to our numerical example, let us consider 
the situation when there is more than one type of coun- 
terion in the solution. Suppose that the solution con- 
tains monovalent (1:1) salt of concentration ci, multiva- 
lent (z:l) salt of concentration Cz, and assume for simplic- 
ity that there is only one species of monovalent co-ions. 
The salt concentrations ci and Cz then determine the 
distribution of ions around the PE - such that far away 
from the PE concentrations of monovalent counterions, 
multivalent counterions and co-ions decay to ci, Cz and 
ci + zcz, respectively. 

Figure 1 b shows results for monovalent salt of concen- 
tration ci — 10 mM and the same 4:1 salt concentrations 
as in part a (solid lines). For each value of Cz there is 
a distinct number of condensed multivalent ions, which 
depends on the multivalent salt concentration. This num- 
ber is very different from the Manning prediction, Eq. J^l, 



for a single counterion species (dashed line). The number 
of condensed multivalent ions is now determined not only 
by a balance of entropy and electrostatics, but also from 
the competition with the monovalent ions. Furthermore, 
this competition can be influenced by ion-specific short 
range interactions leading to a strong influence on the 
number of condensed ions. As a simple example, the dot- 
ted lines show results for multivalent ions that are slightly 
larger than the monovalent ions, having a radius of clos- 
est approach to the PE that is larger by 2 A from that of 
the monovalent ones, leading to a considerable decrease 
in the condensation of multivalent ions. Compared to 
the case of identical ion sizes, the 4:1 salt concentration 
has to be increased by roughly an order of magnitude in 
order to have the same number of condensed multivalent 
ions. 



3. COUNTERION COMPETITION IN DNA 
AGGREGATION 

In Ref. DNA aggregation, induced by multivalent 
ions, was studied experimentally. The conditions for ag- 
gregation were mapped with varying concentrations of 
monovalent salt, multivalent salt, and DNA. These ex- 
periments provide interesting evidence for the role played 
by competition between different ion species. We will 
concentrate on experiments that were done in solutions 
containing the following ingredients: short (150 base 
pair) DNA chains of concentration cdna (base pairs per 
unit volume), spermine (4-|-) of concentration Cz and 
monovalent salt of concentration ci. At sufficiently high 
concentrations the multivalent salt mediates an attrac- 
tive interaction between DNA chains, leading to aggrega- 
tion of the DNA and its precipitation from the solution. 
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ci [mM] 


cl[mM] 


apl 


ap, (PB) 


apz (SR) 


2 


± 0.0003 


0.194 ±0.020 


0.186 ±0.005 


0.191 ±0.006 


13 


0.011 ±0.002 


0.191 ±0.013 


0.178 ±0.002 


0.172 ±0.003 


23 


0.031 ± 0.005 


0.173 ±0.025 


0.172 ±0.002 


0.163 ±0.003 


88 


0.52 ±0.05 


0.135 ±0.026 


0.164 ± 0.002 


0.149 ±0.003 



TABLE I: Excess of 4-valent ions near DNA extracted from DNA aggregation experiments, compared with calculated values 
using Poisson-Boltzmann theory (PB) (treating the DNA as a cylinder of radius 10 A) and Poisson-Boltzmann theory with 
short range interactions (SR). 



For constant ci and with increase of Cz DNA starts to 
precipitate when Cz crosses a threshold value which we 
denote as c^^aggr- The dependence of this threshold on ci 
and cdna was measured in ^ in great detail. 

At the onset of aggregation all the DNA is still sol- 
ubilized. The partitioning of multivalent ions into free 
and condensed ions is thus governed by their distribu- 
tion around isolated DNA chains. Within the relevant 
experimental parameters, ion density profiles associated 
with different chains are decoupled from each other due 
to screening by salt. As a result, a linear dependence of 
Cz, aggr on Cdna is expected theoretically 0: 

C2,aggr = C* + ap*CDNA- (2) 

The first term on the right hand side, c*, is the con- 
centration of multivalent ions far away from the DNA 
chains. This quantity plays the role of the salt concen- 
tration, which determines the distribution of ions near 
the DNA chains. The second term is the contribution of 
condensed ions to the volume-averaged concentration of 
multivalent ions. The coefhcient multiplying cdna, op*, 
is the excess of multivalent ions per DNA base, where a = 
1.7 A is the monomer length (le per a on the chain) and 
pI is a spacial integral of the local excess of multivalent 
ions: 

pl^2^ j dTr[cz{r)-c*z] (3) 

where Cz{r) is the local concentration. Wc will assume 
that c* and p* do not depend on cdna, which is justified 
for sufficiently longchains for which translational entropy 
can be neglected 13- Indeed, the experimental measure- 
ments of Cz.aggr, as fuuctiou of Cdna, fall on straight lines 
(within experimental error bars) across several orders of 
magnitude of DNA and spermine concentrations. The 
coefficients c* and ap* can be extracted from the exper- 
imental data and are reproduced in Table L 

3.1. Comparison with PB theory 

Table I lists the concentrations of monovalent salt, 
ci, multivalent salt c* (extracted from the linear fit to 
Eq. (|21)), and the excess of multivalent ions ap*. The 
last quantity is a measure for the number of condensed 
multivalent ions at the onset of aggregation. Note that 



it is of the same order of magnitude for all four monova- 
lent salt concentrations. However, c* varies by at least 
four orders of magnitude. The large variation in c* is 
a consequence of competition between monovalent and 
multivalent counterions: a relatively small increase in 
monovalent salt concentration requires a large addition 
of multivalent salt in order to keep the number of con- 
densed multivalent ions constant. This point is further 
discussed in Sec. IV. 

Table I provides simultaneous measurements of the salt 
concentrations (ci, c*) and the excess of condensed multi- 
valent ions ap*. Note that no particular model specifying 
the relation between ci, c* and ap* is assumed in Eq. ^ 
and the latter two quantities are obtained independently 
from the linear fit. This data can be used to test any par- 
ticular theory used to calculate ion distributions around 
DNA. Such a comparison, with Poisson-Boltzmann (PB) 
theory, is shown in the fourth column of Table I: for each 
pair of salt concentrations (ci, c*) the excess apz, as cal- 
culated from PB theory, is compared with the experi- 
mental value of ap*. 

For the three smaller values of ci = 2, 13, and 23 mM 
there is a reasonable agreement with experiment (within 
the error bars). However, for ci = 88 mM there is a 30% 
deviation. The overall agreement with PB theory is sur- 
prisingly good, considering that PB theory does not work 
so well for bulky multivalent ions. Ion-ion correlations 
that are ignored in PB theory and are important with 
multivalent ions, tend to increase the number of bound 
multivalent counterions. Instead, for ci = 88 mM, the 
number of bound multivalent ions is decreased. We con- 
clude that ion correlations by themselves are not the main 
source for deviations seen in Table I, and short range in- 
teractions also play a prominent role. 

There are many types of short-range interactions that 
are not taken into account in PB theory. Spermine is 
a long, relatively narrow molecule, which can approach 
DNA at close proximity, and even penetrate the grooves 
at certain sites and orientations Q. On the other hand, 
configurations that are close enough to the DNA are ac- 
companied by a loss of orientational entropy. Other fac- 
tors that modify the interaction of spermine with DNA, 
compared to simplified electrostatic models, include dis- 
persion interactions, specific ordering of charges on the 
spermine and DNA, and arrangement of the surrounding 
water molecules. 

Taking all the above parameters into account is beyond 
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the scope of this work. Instead we demonstrate, within 
the framework of PB theory, that short-range interac- 
tions can influence the competition between monovalent 
and muhivalent ions, and thereby affect the onset of ag- 
gregation in a similar way to that seen it Table I. As 
a simple example (with somewhat arbitrary parameters 
chosen to demonstrate our point) two short-range effects 
are added to the PB model. We consider 4-valent ions 
that are larger than the monovalent ones. Hence the dis- 
tance of closest approach to the DNA is different for the 
two species. In this example these distances are taken 
as 9 A for the monovalent counterions and 12 A for the 
multivalent ones. In addition, we include a short-range 
attraction between the multivalent ions and DNA: multi- 
valent ions gain 3 ksT if their distance from the DNA is 
smaller than 15 A. Qualitatively these are two compet- 
ing effects. The first one (closer approach of monovalent 
ions) slows down replacement of monovalent ions by mul- 
tivalent ions, while the second (short-range attraction) 
has the opposite effect. The balance between the two 
effects is different for different ci and c*. 

The last column of Table I shows values of apz cal- 
culated using the above modified model. These values 
(SR) are shown next to the results of the usual Poisson- 
Boltzmann theory (PB) and compared with the exper- 
imental value of ap*. For ci = 2mM, is almost 
the same in the two calculations. For ci = 88 mM and 
c* = 0.52 mM, pz is considerably decreased with the in- 
clusion of short-range interactions, and is closer to the ex- 
perimental value. Any one of the two short-range effects, 
by itself, results in a large discrepancy with experimental 
data at low salt concentrations. 

We believe that the importance of competing mech- 
anisms for a long, multivalent ion such as spermine go 
beyond the simple modifications to PB described above. 
More refined modifications include the loss of orienta- 
tional entropy at close proximity to the DNA. This effect 
creates a short-range repulsion, whereas the correlation 
effect beyond mean-field is similar to a short-range at- 
traction. Similar competing mechanisms were found in 
simulation of spermidine (3"*") and NaCl in contact with 
DNA 9]. In particular, for high salt concentrations sper- 
midine binding was considerably reduced compared to 
Poisson-Boltzmann theory. In the computer simulation 
m both molecular-specific interactions, the geometrical 
shape of the constituents and ion-ion correlations were 
taken into account. All these effects, and especially the 
geometry of spermidine, which is similar to that of sper- 
mine, were found to play an important role. 



4. TWO-PHASE MODEL FOR COMPETING 
SPECIES 

Two-phase models have been widely used to describe 
the distribution of counterions around cylindrical macro- 
molecules |2|,il0j. In these models ions are considered as 
either condensed or free. The condensed ions gain elec- 



trostatic energy due to their proximity to the negatively 
charged chain but lose entropy, since they are bound at a 
small cylindrical shell around it. For systems with more 
than one type of counterion Manning introduced the so- 
called two- variable theory (l3|, which is an extension of 
his previous model 0, Q . This model has been used to 
analyze condensation (single molecule collapse) of DNA 
molecules induced by spermine and spermidine |12L |13j . 
In this section we present a similar model, which dif- 
fers from Manning's two- variable theory in some details. 
Our main purpose is to explain the large sensitivity of 
c* to changes in monovalent salt concentration. As a by- 
product of our analysis we compare our two-phase model 
with PB theory and Manning's two-variable theory. 



4.1. Model details and main equations 

Assume that the PE is confined within a finite cylin- 
drical cell of area A. The free energy is then written as 
follows: 

+ ^(-PDNA + -ZPz + Pl)0 (4) 

The first two terms are the entropy of bound multivalent 
and monovalent counterions, where pz and pi are the 
number of condensed ions per unit length of the PE. We 
assume that condensed ions are bound on a cylindrical 
shell around the chain and take its area, for simplicity, 
to be 

Ac = n(f (5) 

where d is the radius of the PE. The length A is included 
in order to have a dimensionless argument inside the log- 
arithms, and can be chosen arbitrarily. 

The next two terms are the entropy of free counterions. 
The numbers per unit area of free multivalent ions, pi, 
and of free monovalent ions, p\ , are related to the number 
of condensed ions since the total number of ions within 
the cell is fixed: 

pi = {A~Ac)4^Acz^Pz, 

p{ = [A- Ac)c{ = Aci+ p-DNA- Pi- (6) 

where we introduced the concentrations of free counteri- 
ons c{ and c{. 

Finally, the electrostatic energy is evaluated as if all 
the bound ions are exactly at the cylinder rim, r — d, 
and the linearized Debye-Hiickel approximation is used 
for the electrostatic potential at r > d. This leads to the 
last term in Eq. I^J , where 4> is the electrostatic potential 
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at r = d, given by 

(j) = ~21b{pdna - zpz - pi) 



KdKi{Kd) 



(7) 



where we assume that the outer cell radius is much larger 
than d and , Pdna is the number of unit charges per 
unit length of DNA, Kq and Ki are zeroth and first order 
modified Bessel functions of the first kind, respectively, 
and is the Debye length: 



2c( + z{z + 1)( 



(8) 



The number of condensed monovalent and z-valent 
counterions is found by minimizing the free energy with 
respect to pi and pz, yielding: 



log 



Pz 



clA, 



-zcf) ; log 



Pi 



(9) 



4.2. Consequences for DNA aggregation 

We are interested in the qualitative dependence of c* 
on ci. Note that c* has the same role as c{ in the two- 
phase model. For the monovalent salt, assuming that 
ci > CDNA, Ci can be replaced by ci 0- Equation © 
then yields the following relation: 



£1 
A, 



ciA 
Pi 



(10) 



where p| is the linear density of bound monovalent ions at 
the onset of aggregation . Qualitatively, pi is the only 
ingredient that needs to be estimated in this equation, 
since ci is controlled experimentally and zpl is of order 
one. 

The main outcome of Eq. H1U|I is that c* scales roughly 
as (ci)^. This explains the large variation of c* at dif- 
ferent monovalent salt concentrations since z = 4. There 
are several sources for corrections to this scaling. The 
first one is the dependence of pi on ci and p*. A second 
source of corrections is the effect of short-range interac- 
tions, which was discussed in the previous section within 
PB theory. In addition, Eq. H10(l involves all the approx- 
imations of the two-phase model. 



4.3. Comparison with other models 

We conclude this section by comparing the predictions 
of the two-phase model with those of PB theory and Man- 
ning's two- variable theory (see also This is in- 
structive due to the wide use of two-phase models in the 
literature. Table II lists the value of c* calculated with 
the two-phase model, using the values of ci and p* of 
Table I. The two-phase model can be seen to agree quali- 
tatively with PB theory. Quantitatively, their predictions 
differ by a factor of up to four in the table. 



Cl 




c* (two-phase) 


cl (PB) 




c* (Manning) 


2 


0.194 


4.1 X lO"'' 


4.3 X 10" 


4 


7.5 X lO"'' 


13 


0.191 


1.0 X 10"^ 


3.7 X 10" 


2 


3.4 X 10"" 


23 


0.173 


7.4 X 10"^ 


3.3 X 10" 


2 


4.1 X 10"'' 


88 


0.135 


3.9 X 10"^ 


1.1 X 10" 


1 


4.8 X 10"^ 



TABLE 11: Comparison of two-phase models (two-phase, 
Manning) with PB theory. All the concentrations are in mM. 



Our two phase model differs from Manning's two vari- 
able theory in some details. First, the area used in the 
expression for the entropy of bound counterions is differ- 
ent. Second, the expression for the electrostatic energy 
of bound ions is given in Manning's theory by: 



21b {p 



DNA 



zpz~pi)\og{l-e-^'') 



(11) 



Note that for small nd the two forms in Eqs. Q and 1)11(1 
are similar if d is replaced by a, since: 



j^o(^d) 
ndKi (nd) 



log 1 



^ — K,d\ 



(12) 



In the last column of Table II we present the results 
of Manning's two variable theory, in the version that was 
used in Refs. 

Ill 111 (with different areas of con- 
densation for monovalent and multivalent counterions). 
Compared to our two phase model, deviations from PB 
theory are larger, typically of approximately two orders of 
magnitude. Since both two-phase models are quite simi- 
lar to each other, their different predictions demonstrate 
the large sensitivity to model-dependent parameters. In 
our opinion such models are useful for obtaining qualita- 
tive predictions, but should be used with great care when 
quantitative predictions are required. 



5. SUMMARY 

In this paper we discussed competition between ions of 
different valency in DNA aggregation, concentrating on 
DNA-counterion complexes in the dilute (non-aggregate) 
phase. Due to competition, the number of condensed 
multivalent ions is highly sensitive to salt concentration 
and to short-range, ion-specific effects. Simplified models 
that include only electrostatic interactions are thus lim- 
ited in their capability to predict the conditions required 
for aggregation. An important experimental evidence for 
the importance of specific interactions is that different 
multivalent ions vary strongly in their ability to induce 
condensation or aggregation of DNA, even if they have 
the same valency 0, In addition to the role of 

specific interactions in the dilute DNA phase, they also 
play a prominent role in the aggregates |20l [2ll | , where 
the gap between neighboring DNA chains is typically of 
order 10 A [H. 
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